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We study the macroscopic entanglement properties of a low dimensional quantum spin system by 
investigating its magnetic properties at low temperatures and high magnetic fields. The tempera- 
ture and magnetic field dependence of entanglement from the susceptibility and magnetization data, 
comparing the experimental extraction with theoretical estimates are given. Extraction of entan- 
glement has been made possible through the macroscopic witness operator, magnetic susceptibility. 
The protocol followed in doing so has been outlined in some recent work. Various plots of entan- 
glement comparing the experimental extraction with theoretical estimates are given. Quantitative 
comparison between concurrence and entanglement witness is given for both the theoretical and 
experimental results. Theory and experiments match over a wide range of temperature and field. 
The spin system studied is a chain, which exhibits dimerisation and yields fascinating entanglement 
properties when the temperature and magnetic field are varied. These spin systems exhibit quantum 
phase transition (QPT) at low temperatures, when the magnetic field is swept through a critical 
value. We show explicitly for the first time, using tools used in quantum information processing 
(QIP), that quantum phase transition (QPT) can be captured experimentally using canonically con- 
jugate observables. Macroscopically, quantum complementarity relation clearly delineates entangled 
states from separable states across the QPT. We have estimated the partial information sharing in 
this system from our magnetization and susceptibility data. The complementarity relation has been 
experimentally verified to hold in this system. 
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I. INTRODUCTION 



Quantum comp utat ion and information has had tremendous appeal to researchers, both on the theoretical and 
experimental fronts^l. The most important element for quantum protocols is the quantum correlation arising from 
entanglementP^. Entanglement in quantum condensed matter systems has been extensively studied from the perspec- 
tive of Quantum Information Processing for the last decade. Various measures of entanglement, and its quantification 
are now available^. Recently, there has been significant development in quantification of entanglement, especially 
in spin systems, where extraction of entanglement has been made possible through macroscopic witness operators 
like magnetic susceptibility^^. Spin chains have been proposed for possible applications in making components like 
quantum wires which can be used to connect the various quantum registers or gates in a quantum computer. The 
study of entanglement c haract eristics of spin chains is important for their application in implementation of quantum 
communication protocol d 11 ! 12 ! as well as for their potential application in constructing hardware components of a 
quantum computer. This has led to study of different properties of Heisenberg spin chains that would be used for 
computations as well as quantum communications^. Spin systems, with their finite dimensional Hilbert spaces are 
ideal for studying the dynamics of entangled states. For this purpose the Heisenberg spin Hamiltonian with its ex- 
change interaction is particularly instructive^^. Since for the purpose of teleportation and other quantum protocols, 
maximally entangled channels provide quantum information transfer with ver y hig h fidelity, it is important that one 
should be able to quantify the entanglement of a particular system or channet^M 

In the present work, we have studied the magnetic properties of spin chains, which exhibits dimerisation and yields 
fascinating entanglement properties when the temperature and magnetic field dependence of the magnetic suscep- 
tibility is studied. We have extracted the temperature and magnetic field dependence of entanglement from the 
susceptibility and magnetization data. These spin systems exhibit quantum phase transition (QPT) at low temper- 
atures, when the magnetic field is swept through a critical value. We have also obtained signatures of QPT through 
our witness operators. Here we have estimated the information sharing from our magnetization and susceptibility 
data, which are known to be complementary observables. 

Specifically, we study the entanglement in copper nitrate (Cu(NO^)2 x 2.hH 2 0) using magnetic susceptibility as 
marcoscopic entanglement witness, where we report the variation of entanglement as a function of temperature as 
well as magnetic field. This system exhibits dimerisation and presents us the opportunity to study quantum phase 
transition (QPT) owing to a reduction in the relevant degrees of freedom, enabling one to capture the dynamics in 
a reduced Hilbert space. Due to this reduction, there is a clear crossover of the ground state and the first excited 
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state as a function of magnetic field at low temperatures. We capture essential features of quantum phase transition 
through suitable quantum information processing tools. In this system the Cu 2+ has an unpaired, localized electron. 
It is known to be an alternating dimer spin chain, where the interaction between alternating spins (say spin number 
2 and 3) is much less than that between the other two spins (spin number 1 and 2). This paper has been divided 
into two major sections, one devoted to study of entanglement in this system and the other is devoted to the study 
of quantum phase transition using the tools of quantum information processing. 

In the next section, we start with reviewing the entanglement property of spin chains, especially pertaining to its 
quantification through macroscopic entanglement witness in the form of susceptibility 7 '. Subsequently we proceed to 
the magnetic properties of dimerised spin system, for explicitly demonstrating the effect of temperature and magnetic 
field on entanglement. In the following section, we present the high field magnetization data, especially pertaining 
to the observation of QPT. It is observed in the present case that a QPT separates the entangled and non-entangled 
Hilbert spaces. 
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FIG. 1: (a) Susceptibility data (circles) of Cu(NC>3)2 x and fit (solid line) to Eq. [3] corresponding to the alternating 

dimer model (Eq. [3l. The dotted line marks out the separable regime (Eq. Ill from the entanj;" 
Extraction of entanglement (concurrence) from the susceptibility data of copper nitrate using Eq 
theoretical fit to this experimentally extracted data using Eq. [9] 



111 from the entangled regime (see text), (b) 
flOl The solid line is the 



II. ENTANGLEMENT IN SPIN CHAINS 



We briefly review the entanglement properties of different quantum states of a spin chain in this section. The 
Hamiltonian of a spin chain with nearest neighbour interaction, is given by, 

H = jJ^Si-Si+i + B^Sf (1) 

i i 

where, the spins are given by (in the unit of h—1), Si = = 1,2, 3 are the three Pauli spin matrices mentioned 
earlier. 

The magnetic field is assumed to be applied along the z direction. The net magnetization would be defined as the 
expectation value of 03, the eigenvectors of which are | t) = f 1 an< ^ ^ ^ = 

Since most of the entanglement measures currently available are for bipartite systems, it is worth exploring how 
one can extract the bipartite entanglement from a given susceptibility data for a spin chain. We consider a two spin 
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system. The Hamiltonian for this two qubit system can be written as: 



n = js 1 -s 2 ^b{s z 1 +s z 2 



D 



(2) 



The eigenvalues and the corresponding eigenvectors of this Hamiltonian are: E\ = 4+B; | ^t) 
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where 



\t/j + ) and \ip~) are the Bell states. 

A ferromagnetic system (J < 0), in the absence of external magnetic field doesn't display any entanglement in the 
ground state. This is because the order parameter (magnetization) commutes with the Hamiltonian and there is no 
spin fluctuation at low temperatures. However, for an antiferromagnetic ground state, i.e., for systems exhibiting 
J > 0, the order parameter, which is the staggered magnetization, does not commute with the Hamiltonian and thus 
there will be spin fluctuation, even at zero temperature. Hence, this system is entangled in absence of any external 
field. In fact, the ground state of an antiferromagnetic interaction ( J > 0) is the maximally entangled state 
for the above system, whereas, the ground state for J < 0, is an equal mixture of the three fold degenerate triplet 
states and is a separable state. 

In absence of external magnetic field the order parameter (staggered magnetization) doesn't commute with the 
Hamiltonian leading to non-trivial dynamics at zero temperature. Thus the order parameter can be used as an entan- 
glement witness, an observable whose expectation value is non-zero for an entangled state and is zero for a separable 
state. Magnetic susceptibility has been shown to be a macroscopic entanglement witness^, where the authors show 
that susceptibility is a measure of entanglement, which is a non-local property of the spins. Its canonically conjugate 
observable, magnetization is a measure of the local spin properties and is therefore an observable complementary 
to the magnetic susceptibility. The magnetic susceptibility ca n be measured experimentally and it is possible to 
extract a measure of entanglement from the experimental dat a 1 19 * 20 -*. One of the first systems which exhibited en- 
tanglement through bulk measurements like magnetic susceptibility or heat capacity is the insulating magnetic salt 
LiHo.Yx^FjsS. 




FIG. 2: Theoretical value of concurrence as a function of magnetic field and temperature for a dimer model. The temperature 
and magnetic field values are in arbitrary units. 

The magnetic susceptibility, x, can be defined as the field derivative of magnetization M (in the applied field 
direction, which here is z) in the limit B — > 0, and the magnetization, M, is defined as the sum of the expectation 
values of the z-component of individual spins, M = X)i=i( <7 2:)- The expression for susceptibility derived from an 
alternative dimer model is given a d 19 l 2 ° l, 
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X = (^)=(g 2 A/kBT)(Mi) B=Q = (g 2 f il/k B T)Y,(S l z S* +1 ) « (g 2 fx%N/k B T)[(l/4) + «S X • &»/3], (3) 



where, AT is the total number of spins per mole, p B is the Bohr magneton and k B is the Boltzmann constant. It was 
derived using the isotropy of space, nearest neighbour interaction and the fact that (Si + S 2 )=0 for an antiferromagnet. 
One can evaluate (S 1 • S 2 ), from the bipartite Hamiltonian (Eq. [2J, for the B = limit, yielding, 

= ft • (4) 

To extract the entanglement, one needs to consider the definition given by Wootters^where it has been shown that 
concurrence is a good measure of entanglement. Concurrence is given by, 



C = max(0, v Ai — y A 2 — y A3 — y A 4 ), (5) 

where Ai > A 2 > A3 > A4, are the eigenvalues of the operator, 

P12 = 02 ® ® cr 2 , (6) 

such that, P12 is the two particle reduced density matrix and the asterisk denotes complex conjugation. It has been 
shown by O'Connor and Wootters^, that for an antiferromagnet, the concurrence is given by, 

C = ^max[0,H-l], (7) 

in absence of magnetic field. Here, U is the free energy, which, by definition, is the expectation value of the Hamilto- 
nian, i.e., U — (H) . Thus, considering the system to be isotropic, the formula for concurrence, given in Eq. [7j can be 
written as, 

C = 2 max [0, (S 1 ■ S 2 ) - (1/4)]. (8) 

Using Eq. [4] we have, 

C = 2 max [0, ^— r ] (9) 

1 + 3e kT 

From the condition of separability of two particle states and Eq. ^ one can obtain the concurrence in terms of 
susceptibility, 



Thus for a separable state, the magnetic susceptibility will satisfy, 

g 2 p%N 1 . , 

v > : L±M , (11) 

X ~ kT 6 V ; 

Entanglement witness, which in spin systems is susceptibility is a more general quantification of entanglement and 
hence could be applicable to a general spin chain, with spin - |, spin-1 or spin-|. However, concurrence is defined 
only for spin - h system^! and this measure of entanglement has a striking match with magnetic susceptibility as a 



measure of entanglement. Here, the expression for concurrence given in Eq. 10 is the entanglement witness which 
will be extracted from experimental data below. Higher spin systems will contain higher dimensional entanglement 
and thus the deviation from entanglement witness is more 2 ^. The system that we have considered for this paper is 
copper nitrate (Cu(N03)2 x 2.5i/ 2 0), an inorganic compound and is a spin - \ system. The sample was procured 
from Sigma Aldrich. Cu{NO^)2 x 2.5H 2 has been shown to be an alternative dimer spin chairP^. The corresponding 
spin Hamiltonian would be of the type, 



5 




FIG. 3: Experimental value of concurrence as a function of magnetic field and temperature for Cu(NOz)<2 x 2.5H2O. The 
magnetic field values are in Tesla and the temperature is in Kelvin. 



H — ~^^{J\^>2i ■ ^2i+l + J2&2i+1 ' 821+2), (12) 

i 

representing pairs of spins that are alternately coupled by strong intra-dimer coupling J\ and weak inter-dimer 
coupling J 2 fa 0.25 J-P^. This system behaves like a chain of dimers or exchange coupled pairs which are themselves 
weakly coupled with each other. 

Entanglement in Cu{NO^)2 x 2.5H2O has been reported earlier by Brukner et ai 19 , who used the susceptibility 
data of Berger et aP^ to extract entanglement. We have measured the temperature and magnetic field dependence 
of susceptibility of Cu(NO^)2 x 2.5H 2 from 2 K to 20 K and from to 4 Tesla. The temperature dependence of 
susceptibility in zero field is shown in fig [lja) . All magnetic data was obtained using Quantum Design Magnetic 
Property Measurement System (MPMS). 

The zero field susceptibility data was fitted to the formula for susceptibility given in Eq. [3] (see Figflja)). The 
solid line shows the fitted curve. The data (open circles) and the fit (solid line) are in good agreement with each 
other. The dotted line shows the entanglement boundary given by Eq. the entangled region is represented by 
the region towards the left of the dotted line. The copper spins exhibit entanglement up to the temperature where 
antiferromagnetic correlations persists, i.e. up to around 5 K. This is slightly more than the ordering temperature, 
which corresponds to the peak in the susceptibility curve. Fig. [ljb), shows the extracted value of entanglement from 
the experimental data (open circles) , which were extracted using the expression for entanglement witness (concurrence) 



given in Eq. 10 The theoretical fit (solid line) to this experimentally extracted data is also shown here. It was 
generated using the expression for concurrence given in Eq. [9] for J=4 K. It can be seen in this figure that the 
concurrence vanishes above the temperature where the antiferromagnetic correlations vanish. Here too the theoretical 
fit (solid line) matches reasonably well with the experimentally extracted values of concurrence generated using the 
entanglement witness (fig. [ljb)). Thus entanglement is, not only an entity which is quantifiable, but it also is 
something that can be measured experimentally. 

The effect of temperature and magnetic field on the entanglement of a Heisenberg anti-ferromagnet, considering 
the dimer model has been dealt theoretically by Arnesen, Bose and Vedraf^, wherein they consider the effect of 
temperature induced magnons on entanglement. In fig. [2] we re-generated the theoretical 3D plot given in fig. 1 of 
Arnesen et. aZP^l, however, unlike theirs, which was plotted in arbitrary units, here the temperature and magnetic field 
has been scaled exactly to match with experimental data. The magnetic field values are in Tesla and the temperature 
is in Kelvin. The entanglement (concurrence) in this plot was generated for the dimer model, considering J = 4 K. The 
experimental 3D plot of concurrence as a function of field and temperature is shown in fig. [3j The visual difference 
between the theoretical and experimental 3D plots is due to the fact that the lowest temperature in the theoretical 
plot is zero Kelvin, and the maximum magnetic field is 7 Tesla, whereas the minimum temperature displayed in the 
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FIG. 4: Experimental value of concurrence as a function of magnetic field at T = 2K (open circles). The theoretical concurrence 
at T = 2K is also plotted as a function of field (dashed lines) . 



experimental 3D plot is 1.8 K and the maximum field is 4 Tesla. Though we have collected magnetization data up 
to 7 Tesla, we have shown the experimental plot only up to 4 Tesla. This is partly because most of the relevant 
physics pertaining to this section is captured up to 4 Tesla, and partly because of some subtleties involved in the 
analysis of susceptibility data above 4 Tesla, which is mentioned in the next section. One can see the general pattern 
of decreasing entanglement with increasing temperature is consistent in both the theoretical and experimental plots, 
so is the decrease in entanglement with increase in magnetic field at the lowest measured temperature. This is due 
to the greater contribution from the separable singlet state in the statistical mixture as the temperature is increased 
at a particular field or as the magnetic field is increased at a particular temperature (see fig. 3| . The plot also shows 
that the entanglement vanishes at the same temperature (around T = 5K) irrespective of the strength of magnetic 
field. 

In fig. [4j we have shown the plot of experimental value of entanglement as a function of magnetic field at T=2 K. 
Also shown in this figure is the theoretical entanglement as a function of field. As expected, this shows a gradual 
decrease from the maximum value. To understand this let us consider the following. At a finite temperature, say 
T = 2K, the state is in a mixture of the four eigenstates, the thermal density matrix of this state is given as, 



P = -{l-TXrie^ + + | ttXtt |e-tf- B >" + | UHU \e~^ +B ^}, 



(13) 



where, = ^ and Z = Tr(p) = e^l 3 + + e ^+ B ^ + e ^~ B ^. 

At finite (non-zero) temperatures the system is in a statistical mixture of all the eigenstates. At temperatures as 
low as 2K two lowest energy states, and | 44-), will have the dominant contribution. Thus, as the magnetic field 
is increased, the contribution from | 4-4-): a separable state, increases gradually at the expense of contribution from 
\^p~), which is a maximally entangled state. As one keeps on increasing the magnetic field, the | 4-4-) state ought 
to become cross the \ip~) state and become the lowest energy state above a critical field. This corresponds to the 
quantum phase transition (QPT) as can be seen from the theoretical 3D plot of concurrence ^ as a function of 
field and temperature. Till the quantum critical point (QCP), the eigenstates do not change, though the eigenvalues 
would. At zero temperature the ground state would remain intact below QCP, but for finite temperature we would 
have a statistical mixture, skewed in favour of the two states mentioned above, closing the excitation gap between the 
two states. However, despite the closing of this gap, one can still see signatures of QPT at finite temperatures, and 
the detailed analysis is given in the next section. 

In fig. [4j we can see that the experimental variation of concurrence at a fixed temperature is in reasonably good 
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FIG. 5: (a) Energy eigenvalues vs magnetic field for a dimer system, (b) Energy eigenvalues vs magnetic field for an alternating 
dimcr system. 



agreement with the theoretical curves especially at low magnetic fields. The slight difference between the theoretical 
and experimental curves at high fields (above 3 Tesla) is probably due to slight difference in J in the two scenarios. 
This discrepancy arises because the theoretical model considered by Arnesen et is a dimer model (with J — 4 K) 
whereas Cu(NO^) 
to that, 



)2 x 2.5H2O is an alternative dimer system corresponding to Eq. 12 with J± 
the isotropy of the system i.e the fact that Xx 



0.25 J2. In addition 



Xy = Xz breaks down in the case of high magnetic field. 
However, since this is a dimerised system, leading to formation of singlet states, where the net S = 0, there is perfect 
SU(2) symmetry, as the system is described by isotropic Heisenberg model, one can use the formula for entanglement 
witness till a field of 4 Tesla for a system where J « AK. 



III. QUANTUM PHASE TRANSITION 



Quantum phase transition (QPT) in a quantum spin system takes place when a level crossing occurs between the 
ground state and excited states as one varies the magnetic field, such that an excited state becomes the new ground 
stateP. This can happen either in a transverse field Ising model^ or in isotropic Heisenberg model where dimerization 
occurs^. The point where the cross over takes place is called the quantum critical point. Behaviour of entanglement 
at the quantum critical point has been extensively studied theoretically by several groups^HUl They have found that 
the entanglement scales close to the quantum critical point (QCP). This means that the entanglement, which in turn is 
related to the correlation function, is scale invariant close to the QCP. The QPT in dimer like systems is characterised 
by crossing of energy eigenstates as a function of some parameter like the magnetic field. This scenario is depicted in 
fig. |5ja) for a dimer system (two qubit system), where the state | ft) crosses the ground state at some critical field 
and replaces it as the ground state, thus changing the symmetry of the ground state. Thus the ground state crosses 
over from a maximally entangled state to a separable state for a two qubit system. The exchange interaction for this 
system, J/ks — 4 K. This is chosen to suit copper nitrate. 

The minimum amount of spins needed to fully describe an alternating dimer system is four, the dimension of the 
Hilbert space being sixteen. The Hamiltonian for such a 4-qubit system is given by, 



T~L — J\S\ ■ S2 + J2S2 ■ 5*3 + J1S3 ■ S4 + B(Si + S*2 + S3 + S4), 



(14) 



where J\ and J2 are the two exchange interactions, with J\ being significantly larger than J2. The energy eigenvalues 
for this Hamiltonian has been plotted as a function of field for Ji/ks = 4 K and J2 = 0.25 x J\. One can see that 
there are two crossovers in the ground state, as is depicted in fig. [5](b). Thus the system will have a lower critical 
field and an upper critical field. 
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FIG. 6: Experimental data of magnetization collected at T = 2 K for Cu(NOg)2 x 2.5H2O (open circles) as a function of 
magnetic field and the fit to the theoretical curve (solid blue line) derived using the alternating dimer model. The dotted line 
(red) is fit to the experimental data using a dimer model (Eq. 16 1. 
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FIG. 7: Experimental data of magnetic susceptibility collected at T = 2 K for Cu{NOz)2 x 2.5H2O (open circles) as a function 
of magnetic field and the fit to the theoretical curve (solid line) represented by Eq. 
experimental data. 



18 The dotted line is a guide to the 



In a significant effort to understand the quantum correlations close to a QPT, concepts from quantum information 
processing, especially that of entanglement was invoked^. This was done especially to extract the extra correlations 
present in QPT, which do not have a classical counterpart. It was found that close to the critical point, the entangle- 
ment depends strongly on the magnetic field. Now since entanglement measures the non-local effect in the system, 
we also tried to capture the non-locality in the spin system close to QCP. 

Here we present the detailed analysis of the magnetic field dependence of magnetization and susceptibility. It 
has been reported by Wiesniak et. oiJ^ that magnetization and susceptibility are complementary variables. Thus, 
whereas magnetization captures the local properties of spins, susceptibility exhibits non-local features in a quantum 
spin system (as is evident from its relation with correlation function). The non-local nature of entanglement intuitively 
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points towards the connection between entanglement and susceptibility (c.f. Eq. 10 ) which in turn is connected with 
spin-spin correlation functions. 

Wiesniak et. alV^ have derived a macroscopic quantum complementarity relation, which using their notation is 
given as, 



kTx , (M) 5 



Ns y N 2 ^ 

non-local properties local properties 



< 1. (15) 



Again following their notation, let Q=l — and P = , where P describes the local properties of individual 
spins, Q the non-local quantum correlations between spins. Since Q, through its dependence on the susceptibility is 



proportional to two-point correlation function, its positive value implies the existence of entanglement (c.f. Eq. (10)). 
When Q = 1, the system is maximally entangled and in such a case, the local property of spins, P = 0. In the 
other extreme, when P — 1 or the magnetization is a maximum, the local properties of the spins is well defined at 



the expense of the non-local property, i.e., Q = 0. Thus, mathematically, the above relation (151 describes partial 
quantum information sharing between local and non-local properties of spins. Following the work of Wiesniak et. al., 
a similar complementary relation was also demonstrated theoretically by Tribedi et. aZp^for a spin 1 system. 

At the heart of this inequality lies the fact that the spin operators, S x , S y and S z , or the Pauli spin matrices (in 
case of a qubit) do not commute. The inequality, (cT x } 2 +(a y ) 2 +(a z ) 2 < 1, where the average is taken over an arbitrary 
state, suggests that for a mixed state the spin lies inside the Bloch sphere, whereas for a pure state it will lie on the 
Bloch sphere. In the latter case, the equality relation holds as in such a case it is a pure state. For a particle with spin 
s, this inequality reads, (S x ) 2 + (S y ) 2 + (S z ) 2 < s 2 . If we consider (S 2 } 2 + (S 2 ) 2 + (S 2 ) 2 , the value will always be more 

than s 2 and will in fact be s(s + 1). This extra value comes from the fact that S x , S y and S z do not commute. Thus, 
there is an uncertainly in determining the values of all three components of spin simultaneously. For a composite 
system, this results in extra correlations which essentially is non-local. For such a composite system, entanglement 
plays the role analogous to that of a one qubit coherently superposed pure state that lies on the Bloch sphere. Since, 
the susceptibility, x ^ ((M 2 ) - (M) 2 ), the second term being less than or equal to N 2 s 2 , N being the total number 
of spins, whereas the first term being proportional to A 2 s(s + 1), the inequality given in Eq. 6 of Wiesniak et.al^ 
holds. 

We have analysed our magnetization and susceptibility data extensively and tried to ascertain the experimen- 
tal validity of this relation. In what follows is a careful exposition of metrology and the validity of susceptibility 
measurements at high magnetic fields. 

The magnetization data for copper nitrate is given in fig. [6] Though the entire field and temperature range has 
been shown in fig. [3j here we focus on the data set below the antiferromagnetic ordering temperature, where all the 
interesting features are occurring. The data was fitted to both the dimer as well as the alternating dimer model. Eq. 
[16] gives us the expression for a dimer model. 

M-Tr{p-S z ) - - - 2cos h(2gn B Hk B T) + exp(J/k B T) ' ^ 

where "g" is the Lande g factor. 

The fit to this expression is given in blue. We can see that the dimer model is not a very good representation of 
the real spin system. The expression for alternating dimer model was generated numerically and subsequently used 
to fit the magnetization data. The expression is not explicitly shown here owing to its cumbrous nature. We can see 
that it is a very good fit. The value of "J" obtained from this fit to all the four data set shown in fig. |6]is around 4K. 
This only goes to show that the Cw(A r O--j) 2 x 2.5H?0. is an alternating dimer system as is also evident from neutron 
datcP. 

Hereafter, we analyze the susceptibility data as a function of magnetic field acquired at the minimum temperature 



(2 K), shown in fig. 18 The susceptibility data shows a down turn after 4 Tesla . Here it is imperative to mention 
a few lines about the nature of measurement of generalized susceptibility employed in the MPMS or other magnetic 
susceptibility measurement system. 

The zero field magnetic susceptibility is defined as, \z — (dM z /dH z ). Now from the fluctuation dissipation theorem, 
the magnetic susceptibility is given by: 

X, = ^A 2 (Af,) = ^«M|> - (Af,) 2 ) = -,V- | >>>■■!> (>;•<) I- (17) 
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FIG. 8: (a) Plot of theoretical P+Q as a function of magnetic field and temperature using f 53^=1 ( s z s z) — s\j J for 

constructing the susceptibility (b) The same where susceptibility was constructed using [J2i'j=i{ s 'z s i)) ■ 
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FIG. 9: Side view of the theoretical P+Q as a function of magnetic field and temperature shown in fig. |8|b). 

where A 2 (M Z ) is the variance of magnetization. This definition of susceptibility gives a measure of the correlation 
function, (s*s-£), which could have a magnetic field dependence and one can measure the field dependence of 
susceptibility and find the field dependence of correlation function. 

However, in order to measure the generalized magnetic susceptibility, as a function of magnetic field (and not zero 
field susceptibility), the same procedure could be followed, but the value obtained is the derivative of magnetization 
with respect to field. This approximation mentioned in the equation above, breaks down above the quantum critical 
point, where, B > J/ks and the assumption of isotropy of space breaks down. The validity of this assumption 
holds for levels of B which is smaller or comparable with J such that the system is still to a large extent isotropic. In 
addition to the isotropy of space, there are other subtleties involved in the extraction of relevant quantum information 
from the susceptibility data and is indicated below. 
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FIG. 10: Theoretical susceptibility as a function of field at temperature T =0.5 K, using the correlation function 
(Ei^=i<44) - (Eili 4> 2 ) ('+')• Tne ter m (e£Li s lf) at T =°- 5 K is shown in the fi S ure (dotted curve). The curve 
'*' shows the susceptibility as a function of field at T =0.5 K, obtained using ^Efj^i^z 5 ^)) • This is also obtained by adding 
the first two curves. The curves with the symbol 'o' represents 'P' and 'x' represents 'Q' at T =2 K, where the latter uses the 
susceptibility plotted in the curve '*'. 



We have obtained a theoretical relationship for the field dependence of susceptibility and hence the correlation 



function and is given in Eq. 18 



X(H) = ( 5 VI^)(0.75 + 46 ( ., +H)fj I , J)0 I 3J / 4 (J m ), (18) 

wher^/? = 

Fig. [7] shows the fit to the experimental susceptibility data as a function of field, collected at 2 K. Here we have 
shown the plots only up to 4 Tesla. The experimental value of susceptibility as a function of field (open circles) 
matches the theoretical value (bold curve). However, above the quantum critical field, there is a down turn in the 
data (see fig. [7]) as using the experimental technique employed here to measure the susceptibility, one is not measuring 
the magnetic field dependence of the correlation function, X)t=i( s z s z)- I n this limit, these two are not the same, i.e. 
x(H) = j^p (Sfj=i( s l s z) — (J2iLi s z) 2 ) 7^ dM z /dH. The susceptibility measurement measures only the latter and 
this was also verified by taking the derivative of the magnetization data (dM/dH). 

The only way to ensure that we are really measuring x{H) — Ei j=i( s l s i) an d not the derivative dM/dH, is to 
measure the correlation function using inelastic neutron scattering measurements. 

Another important point is that in a dimer like antiferromagnetic system, since the average spin, (s l z ) — 0, at 
least below a critical field, where the ground state is the singlet Bell state \^~)- However, above a critical field 
where the quantum phase transition takes place, the ground state crosses over to a triplet like state, which is a 

product state and (s l z ) ^ 0. In the definition of correlation function, ^X)fj=i( s z s z) — (^2iLi s z^ V the second 

term will start contributing to the susceptibility, after the quantum critical point. Thus, when we use the term 

(Sf^i ( s z s i) ~ (J2iLi s z) 2 ) m the susceptibility to plot the P+Q as a function of field and temperature, it shows 

an asymmetry either side of the quantum critical point even at the lowest temperature. This is shown in fig. 8ja), 
which gives the theoretical 3D plot of P+Q as a function of field and temperature. This is what one would obtain 
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FIG. 11: (a) Plot of experimental "P+Q" as a function of magnetic field at T =2 K. Here the "Q" used was obtained considering 
the corrected susceptibility, (b) Plot of theoretical "P+Q" as a function of magnetic field at T =2 K. Here the "Q" used was 
obtained considering the corrected susceptibility. 



if we consider the experimental susceptibility in constructing Q = 1 — . This asymmetry is in disagreement with 
what Wiesniak et. al^ obtained. If we ignore the second term in constructing susceptibility, i.e. if we consider only 
(j2i ! j=i( s z s i)^Ji then we obtain the 3D "P+Q" plot as obtained by Wiesniak et. al. and is shown in fig. 8m). 

Before constructing the 2D "P+Q" as a function of temperature, one needs to demonstrate the different zD curves 
of susceptibility mentioned above as a function of temperature. Fig. [10] shows the susceptibility as a function of 

field at T =0.5 K, constructed by using the correlation function, fX}ij=i( s z s z) ~ (J2iLi s z) 2 ) (curve "+"). The term 

(j2iLi s l) 2 ) i s a l so plotted (dotted curve) at the same temperature as a function of field. The sum of these two curves 
is also shown in this figure (curve "*"). Next we show the non-local quantity "Q" as a function of magnetic field at 
T =0.5 K (curve "x"). This was constructed using the modified susceptibility (i.e. using the curve "*" mentioned 
above). The local term "P" is also plotted in this figure (curve "o") at the same temperature as a function of magnetic 
field. We have chosen T =0.5 K for all the 2D plots for consistency as we wanted to demonstrate the difference in the 
susceptibilities owing to different terms in the correlation function. This temperature T =0.5 K was chosen to bring 
out the differences clearly as both " P" and " Q" are quite sharp at such a reasonably low temperature. 



Next, we show the theoretical plot of "P+Q" as a function of magnetic field at T =2 K in fig. 11 V). It shows a 
dip at the quantum critical point which is between 3.5 Tesla and 4 Tesla. To match the experimental value of "P+Q" 
with the theoretical value at T =2 K, we consider the fit to the experimental susceptibility shown in fig. [7] Here 
in constructing "Q", we consider the susceptibility data till the magnetic field where the theoretical curve matches 
the experimental data and then extrapolate the data up to 7 Tesla, using the fitted curve. Since the experimental 

data measures the term (^J2iLi s z) 2 ) as weu ; we need to add this term to match the experimental data with the 

theoretical value of Wiesniak et. al.. This plot is shown in fig. [TT|a). The local property "P" was considered using 

the experimental magnetization data at T =2 K using the formulae, P = N2 2 ■ One can clearly see the dip in 



the value of "P+Q" in fig. 11 'a), corresponding to the QCP at around 3.5 Tesla. This suggests that the "P+Q" 
prescription proposed by Wiesniak et. al. indeed is a detector of quantum phase transition. Our experimental results 
also demonstrate the quantum information sharing and complementarity between the two observables, susceptibility 
and magnetization. The "P+Q" inequality given in Eq. [15] is not violated in our system and the minimum value is 
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indeed reached at QPT. 



IV. CONCLUSION 



We have studied a Heisenberg spin chain compound known to exhibit a dimer or alternative dimer like charac- 
teristics. The system has an antiferromagnetic coupling and therefore exhibits an entangled ground state. This 
entanglement which was extracted from the susceptibility data decreases when the temperature is varied from the 
lowest attainable value (1.8 K) and vanishes around 5 K, where the antiferromagnetic correlations die down. We 
also studied the magnetic field dependence of the magnetization and susceptibility and they both fit reasonably well 
to an alternating dimer model. We extracted the magnetic field dependence of entanglement and observed that the 
concurrence decreases with increase in field, owing to an increase in occupation of the excited states in comparison 
to the antiferromagnetic ground state. At the critical field where, there is a crossover between the ground and 
excited states, bringing about a change in symmetry of the ground states, a quantum phase transition occurs. We 
explored the complementarity in the two observables, susceptibility and magnetization, which capture the non-locality 
and local nature of correlations in this spin system. Following an established prescription, the susceptibility and 
magnetization were suitably combined to form operators "P" and "Q", such that the minimum of the sum "P+Q" 
when plotted as a function of magnetic field captures the quantum phase transition. From this measure ("P+Q"), 
we could capture the QPT as a function of field at 2K. This shows that "P+Q" is a good measure for detecting 
quantum phase transitions, even at non-zero temperatures. This could have tremendous implications in quantum 
information processing as entanglement is known to scale close to QPT from theoretical results. This scale invariance 
of entanglement coupled with the "P+Q" measure could be harnessed to execute quantum communications in spins 
chains. In future, we are making an effort to explore the quantum phase transition at lower temperatures. We also 
intend to do specific heat measurements to capture other aspects of entanglement and QPT in this system. 
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